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Abstract 

We present results for heavy-light decay constants, using both propagating quarks 
and the static approximation, in 0(a)-improved, quenched lattice QCD. At /? = 6.2 on a 
24^ x48 lattice we find fo = 185 ^3 (stat) (syst) MeV, fB = l60^l ^ ^ MeV, 

fos/fo = 1-18 ^ 2 ^^"^ Jbs/Ib = 1-22 ^ 3, in good agreement with earlier studies. 

From the static theory we obtain fj^^^ = 253 ^ ^ ^^^^ MeV. We also present results 
from a simulation at (3 = 6.0 on a 16'^ x 48 lattice, which are consistent with those at 
(3 = 6.2. In order to study the effects of improvement, we present a direct comparison 
of the results using both the Wilson and the improved action at /? = 6.0. 



1 Introduction 



The leptonic decay constants, fp, of pseudoscalar mesons composed of a heavy and a hght 
quark play an important role in weak- interaction phenomenology. In particular J'b, or more 
strictly fsV^B (where Bb, the B parameter of mixing, is expected to be close to one), 

is one of the principal unknown quantities needed for the determination of the CP- violating 
phase in the Standard Model, as well as other properties of weak decays. Lattice QCD offers 
the opportunity for a non-perturbative computation of the operator matrix elements which 
are necessary for the determination of the decay constants and B parameters. 

During the last few years there have been several lattice computations of the decay constants 
of "heavy-light" pseudoscalar (and vector) mesons. The results for the decay constant of 
the D meson, obtained using the Wilson action for the quarks, are in the region of 200 MeV 
(using a normalisation for which /^r =132 MeV). For example, in his 1989 review S. Sharpe 
quoted ^ fo — 180 ± 25(stat) ± 30(syst) MeV as his summary of the lattice results. More 
recent simulations with Wilson fermions also give results in this range |@|-|^. The experi- 
mental bound is fn < 290 MeV 

In the heavy-quark limit the scaling law for the decay constant of a heavy-light pseudoscalar 



meson is fp\JMp ~ constant (up to mild logarithmic corrections). Lattice simulations using 
heavy-quark masses in the charm region indicate that there are large corrections to this 
scaling law (of order 40% at the charm quark mass, decreasing to about 15% at the mass of 
the bottom quark) [0, ^ |[. The value of the decay constant of the B meson deduced from 
these simulations is in the region of 180 MeV. The conclusion that there are violations of the 



scaling law is supported by the large value for fp^/Mp deduced from simulations obtained 
using a static (i.e infinitely-massive) heavy quark |^, ||, ^, |T0|, [11], |T^. 

The important results and conclusions quoted above were obtained from simulations in 
which the mass of the heavy quark is large in lattice units (up to about a half). One 
may therefore worry that discretisation errors significantly contaminate the results. In this 
paper we present the results for decay constants of heavy-light mesons computed using the 
0(a)-improved lattice action proposed by Sheikholeslami and Wohlert |jl3|, with which the 
discretisation errors in operator matrix elements (and hence in the computed values of the 
decay constants) can be reduced from 0{mQa) to 0{asTnQa), where mq is the mass of the 
heavy quark [|14| . This formal reduction in discretisation errors provides an important check 
on the stability of results and conclusions obtained with Wilson fermions. 
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The results presented in this paper were obtained from two simulations of quenched QCD, 
using the Sheikholeslami-Wohlert (SW) or "clover" fermion action for the quarks (see Sub- 
section below). Our main results come from a simulation on a X 48 lattice at /? = 6.2, 
for which 60 gauge field configurations were generated. Details of this simulation and the 
determination of the values of the Wilson hopping parameter corresponding to the chiral 
limit, Kcrit, and to the mass of the strange quark have been presented in ref. [|15|. The heavy 



quarks have masses in the region of the charm quark mass and we study the behaviour of the 
decay constants with the mass of the heavy quark. Interpolating to the mass of the charm 
quark itself, and extrapolating the results to the mass of the b quark, we find that our best 
results for the decay constants of the B and D mesons are 

= 185+3 (^*^*) - ^7 (^y^*) (1) 
Ib = 160 + ^ + II MeV (2) 

fPs 110+2 

/d 



1.18 Z 2 (3) 
1.22 ^ J. (4) 
The details of this calculation and a complete set of results are presented in Section |[ 



fB 



The second simulation is on a 16^ x 48 lattice at (3 = 6.0, using 36 configurations. The 
results, which are consistent with those mentioned above, are presented in Section ^. In 
order to study the effects of improvement on the calculation of heavy-light decay constants, 
we have repeated the computation for both the Wilson and SW actions using a subset of 16 
of these configurations. The results and a discussion are presented in Subsection |3^ . 

We have also computed Jb in the static approximation (in which contributions of 0(l/mb) 
are neglected). A discussion of the calculation and of the results is presented in Section ^. 
The result from the simulation at /3 = 6.2, on 20 of the 60 configurations, is 

/r = 253 + II (Stat) + (syst) MeV, (5) 

and the result at /3 = 6.0 on all 36 configurations is 

/r = 286 + + % MeV. (6) 



Finally, Section ^ contains our conclusions. 
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1.1 Improved Action and Operators 



The SW action is 



X,fJ,,U 



XI 



(7) 



where Sp^ is the Wilson action: 



S^ = T.{ li^hi^) -^T. - l,)U,{x)q{x + A) + q{x + + j,)Ul{x)q{x)] I . (8) 

The decay constants of heavy-hght pseudoscalar and vector mesons are computed using 
lattice axial and vector currents as interpolating operators. In order to obtain 0(a)-improved 



matrix elements we use "rotated" operators JT^ 



Q{x){l + ^T 5) r (1-^7- B)g(x), 



(9) 



where F is one of the Dirac matrices (either ^^^^^ or 7^) , and Q and q represent the fields of 
the heavy and light quark respectively. 

In the static effective theory, in which the heavy propagator is expressed in terms of the 
link variables [0, it is sufficient to rotate the light-quark fields only [|1^], i.e. to use the 
operators 

Q{x)r{l-^rD)q{x), (10) 
in order to eliminate the 0(a)-discretisation errors. 



1.2 Renormalisation Constants Zy and Za 

In order to determine the physical values of the decay constants from those obtained in 
lattice simulations using the interpolating operators in eq. (^, it is necessary to know the 
corresponding renormalisation constants. These are defined by requiring that ZaAJ^^" and 
ZyV^^^^ are the correctly normalised currents, where the superscript "latt" denotes that 
the operator is a lattice current. These renormalisation constants have been calculated at 



one-loop order in perturbation theory for the SW action with rotated operators [|T8|: 

Zv = 1-0.10/ (11) 
Za = I- 0.02/. (12) 



In this paper they are evaluated using the "boosted" coupling suggested in ref. JTOl; specifi- 
cally, we use g'^ = 6/(/9mq), where Uq is a measure of the average link variable, for which we 
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take uo = (S/tcrit) ^- It has been suggested EOl that the use of such an effective couphng 



rather than the bare lattice couphng, resums some of the large higher-order corrections, and 
in particular some of the tadpole diagrams. Using the measured values of Kcth from our 
simulations, we obtain Za — 0.97 (0.96) and Zy — 0.83 (0.82) for the simulation at j3 = 6.2 
(6.0). 

In a recent non-perturbative determination of these renormalisation constants, obtained by 
requiring that the correctly-normalised currents obey the continuum Ward Identities, it was 
found that Zy = 0.824(2) and Z^ = 1.09(3) These results were obtained from a 

simulation at /3 = 6.0 for one value of the quark mass. It remains to be checked that the 
results are independent of the quark mass and insensitive to small variations in (3. For 
this reason, we use the perturbative values, given above, throughout the paper. We note. 



however, that the non-perturbative value of Z^ may be larger by about 15%. In ref. [|T5 
we obtained fn/rrip = 0.138 ^ g, using the perturbative value of Z^- A larger value of Za, 
such as Za = 1.09, would bring this result closer to the physical value of 0.172. However, 
we also observed that fx/ f-K-, which does not require Za-, was in very good agreement with 
the experimental value, and therefore we quote values for the ratios fo / fn and fs / fn in the 
following sections. 

The normalisation of the axial current in the static effective theory is discussed in Section 11. 



1.3 Error Estimation 



Statistical errors are obtained from a bootstrap procedure This involves the creation 
of 1000 bootstrap samples from the original set of configurations by randomly selecting 
configurations per sample (with replacement). Correlators are fitted for each bootstrap 
sample by minimising x^- During the fits, correlations among different timeslices are taken 
into account, whereas correlations among different values of the quark mass are neglected. 
The latter correlations are preserved by using the same sequence of bootstrap samples at each 
quark mass. When extrapolating our results to the chiral limit and physical meson masses, 
the correlation matrix for the fitted quantities is estimated from the full bootstrap ensemble. 
All quoted statistical errors are obtained from the central 68% of the corresponding bootstrap 
distributions |2^ . 



We attempt to quantify the systematic error arising from the uncertainty in the value of 
the lattice spacing, a, determined from properties of light hadrons, and from the string ten- 
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sion |T5|. The differences between results obtained using our central value for [GeV] and 
our upper and lower estimates are quoted as systematic uncertainties in the final estimates 
for decay constants in physical units. Hereafter, where we quote two errors, the first is 
statistical and the second is systematic. 



1.4 Extended Interpolating Operators 

In order to isolate the ground state in correlation functions effectively, it is useful to use 
extended (or "smeared" ) interpolating operators for the mesons. In particular, in the static 
theory it has been found to be essential to use smeared operators in order to obtain any 
signal for the ground state p4|. In this study we use gauge- invariant Jacobi smearing on 
the heavy-quark field (described in detail in ref. p5|), in which the smeared field, Q^{x,t), 
is defined by 

QS{x,t)=J2J{x,x')Q{x',t), (13) 

x' 

where 

N 

J(x,x') = E (14) 

n=0 

and 

3 

A(x, x') = ^{5^/,^_£f//(f -i,t) + Ss',x+iUi{x, t)}. (15) 

i=l 

Wuppertal smearing [^, which uses the operator (1 — k^A)^^ as the kernel of the smearing, 
corresponds to = oo, provided that ks is sufficiently small to guarantee convergence. 
Following the discussion in ref. p5[, we choose ks = 0.25 and use the parameter N to 



control the smearing radius, defined by 

2_ S.-|x|V(x,o)p 
Ex\J{x,o)\^ ■ 



(16) 



The values of and r used in each of the calculations below will be quoted in the corre- 
sponding sections. 



2 Decay Constants from the Simulation at p = 6.2 

In this section we present the results obtained for the decay constants of heavy-light mesons 
from our simulation on 60 configurations of a 24^ x 48 at /3 = 6.2, using the SW action in the 
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quenched approximation. The computations are performed for four different values of the 
mass of the heavy quark, corresponding to Kh = 0.121, 0.125, 0.129 and 0.133, and for three 
values of the mass of the light quark, corresponding to ni = 0.14144, 0.14226 and 0.14262. 
The mass of the charm quark corresponds approximately to Hh = 0.129. The value of the 
hopping parameter corresponding to the mass of the strange quark is Ks = 0.1419 + \ and 



the critical value is Kent = 0.14315 + 1 fig. 

The decay constants are determined by computing two-point correlation functions of the 
form 

c^SW = E(o|^?(^)4''(o)|o) (17) 

X 

where Ji and ^''^^ interpolating operators which can annihilate or create the pseudoscalar 
or vector meson being studied. The labels Q, R denote whether a local (L) or smeared [S) 
interpolating operator is being used. In this simulation we use Jacobi smearing with N = 75, 
corresponding to a smearing radius of r = 5.2. The decay constants are obtained from the 
matrix elements of the local operators, which are determined by computing both the C^^ 
and C^^ correlation functions. 

In order to determine the decay constant, it is necessary to know the value of the lattice 
spacing in physical units. This can be done by relating the lattice measurements of some 
dimensionful quantity to its physical value, e.g. the mass of a light hadron or Among 
the other frequently-used choices are the string tension, a/k", and the IP — IS mass split- 
ting in charmonium. Using nip to set the scale in our study of light hadrons [|l5l we found 



a~^{mp) = 2.7(1) GeV, and a mass spectrum in physical units which was close to experimen- 
tal values. Furthermore, our determination of the string tension gave = 2.73(5) GeV. 
Encouraged by the consistency of these results, we use 

a-^ = 2.7GeV. (18) 



However, the study described in ref. |]T3| showed that the measurement of the pion decay 
constant gave a higher value for the scale, i.e. a'^i^f-,,) = 3.4 ^ ^ GeV, using the perturbative 
value for Za- In order to get an estimate of the systematic uncertainties in the final numbers, 
we evaluate all our results using the central value of a~^{fT,) as well, and quote the difference 
as the upper systematic error on decay constants. The lower systematic error is obtained 
from the uncertainty of —0.1 GeV in a~^{mp). 

In an attempt to reduce the systematic errors associated with the value of the renormalisation 
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constant of the axial current, we also compute fD,B/ fn, and determine fD,B by using the 
physical value of /j^. 



2.1 Decay Constants of Pseudoscalar Mesons 

In order to determine the pseudoscalar decay constants, we start by fitting the two-point 
correlation function 

cl>Ut) = E(o|^^(^,t)pt^(o)|o) 

X 

^ ^exp(-MpV2)cosh(Mp(V2-t)), (19) 

where P is the pseudoscalar density, Zps = {0\P^{0)\P) and Lf is the temporal extent 
of the lattice. This correlation function gives the best determination of the masses of the 
heavy-light pseudoscalars. Symmetrizing in Euclidean time, the fitting range was chosen to 
be 13 < t < 22 for all three values of the light-quark mass. Good plateaus in the effective 
mass were observed and stable fits obtained. The values of the masses of the pseudoscalar 
and vector mesons for the twelve nt-i^i combinations are presented in Table |l|. The values 
obtained by linear extrapolation to the chiral limit for the light quark are also tabulated. 

At large values of t, the ratio of correlation functions 



cmt) (oip^(o)iP) 



tanh(Mp(Lt/2-t)) (20) 



is used to extract the pseudoscalar decay constant, where A4 is the temporal component of 
the axial current. The ratio is fitted in the range 15 < t < 22 with the pseudoscalar mass Mp 
(in each bootstrap sample) constrained to its value extracted from fits to eq. (plQ]). In Fig. ^ 
we plot the ratio of correlators together with the fit to eq. ([20| ) as a function of t. Using the 
value of Zps obtained from the fits of eq. (|T9|), the matrix elements of the local axial current 
are obtained. Although there are other ways of determining these matrix elements, we find 
that the ratios in eq. (^) give the most precise determination. 

In Table ^ we present the results for the decay constants (in lattice units) of the pseudoscalar 
mesons for the twelve Hh-i^i combinations, as well as the values obtained by linearly extrap- 
olating the results to the chiral limit. We also tabulate the results for the quantity jp\fMp 
which, in the heavy-quark limit, is independent of the mass of the heavy quark (except for 
a mild logarithmic dependence). 
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Mp 


Mv 




0.14144 


0.924 


+ 3 

— 1 


0.944 


+ 4 

— z 


0.121 


14226 


0.900 


+ 3 
— 2 


0.920 


+ 4 
— 3 






0.890 


+ 4 
_ 3 


0.909 


+ 5 
— 4 




^crit 


0.875 


+ 4 
_ 3 


0.894 


+ 6 
— 4 




0.14144 


0.822 


+ 3 

— 1 


0.847 


+ 4 

— z 


0.125 


14226 


0.799 


+ 3 
_ 2 


0.823 


+ 4 
_ 3 




0.14262 


0.789 


+ 4 
_ 2 


0.811 


+ 5 
— 4 




^crit 


0.773 


+ 5 
— 2 


0.797 


+ 6 

— 4 




0.14144 


0.715 


+ 3 

— 1 


0.745 


+ 4 
_ 2 


0.129 


0.14226 


0.691 


+ 3 
- 2 


0.721 


+ 4 
- 3 




0.14262 


0.681 


+ 4 
- 2 


0.711 


+ 5 
- 4 






0.665 


+ 5 
- 2 


0.695 


+ 6 

- 4 




0.14144 


0.599 


+ 3 
- 1 


0.637 


+ 3 
- 2 


0.133 


0.14226 


0.574 


+ 3 
- 2 


0.613 


+ 4 
- 3 




0.14262 


0.564 


+ 4 
- 2 


0.603 


+ 5 
- 4 




'^crit 


0.546 


+ 4 
- 2 


0.588 


+ 6 
- 4 



Table 1: Masses (in lattice units) of the pseudoscalar and vector mesons 
for the twelve Hh-f^i combinations aX (3 = 6.2 on a 24^ x 48 lattice. Also 
presented are the values obtained by linear extrapolation to the chiral limit 
{ki ^ fi:crit = 0.14315). 
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fp/ZA 




1 l( f,. 






0.14144 


0.086 + 


1 

1 
i 


0.083 


+ 2 

— 1 


0.124 


+ 2 

- 3 


0.121 


14226 


0.079 + 


1 
1 


0.075 


+ 2 
_ 1 


0.116 


+ 3 
- 3 




0.14262 


0.076 + 


2 
2 


0.071 


+ 2 
_ 2 


0.111 


+ 3 
- 3 




^crit 


0.071 + 


2 
1 


0.066 


+ 2 
_ 1 


0.105 


+ 4 
- 4 




0.14144 


0.086 + 


1 

1 
i 


0.078 


+ 1 

— 1 


0.141 


+ 3 
- 3 


0.125 


14226 


0.079 + 


1 
1 


0.070 


+ 1 

_ 1 


0.132 


+ 3 
- 3 




0.14262 


0.076 + 


2 
2 


0.067 


+ 2 
_ 2 


0.128 


+ 4 
- 4 




^crit 


0.071 + 


2 
1 


0.062 


+ 2 
_ 1 


0.123 


+ 4 
- 4 




0.14144 


0.085 + 


1 
\ 


0.071 


+ 1 

^ 


0.163 


+ 3 
- 3 


0.129 


0.14226 


0.078 + 


1 
1 


0.065 


+ 1 
- 1 


0.155 


+ 3 
- 3 




0.14262 


0.075 + 


2 
1 


0.062 


+ 1 
- 1 


0.151 


+ 4 
- 4 




'^crit 


0.071 + 


2 
1 


0.057 


+ 2 
- 1 


0.146 


+ 5 
- 5 




0.14144 


0.082 + 


1 
1 


0.063 


+ 1 
- 1 


0.193 


+ 3 
- 4 


0.133 


0.14226 


0.076 + 


1 
1 


0.057 


+ 1 
- 1 


0.186 


+ 3 
- 4 




0.14262 


0.073 + 


1 
1 


0.055 


+ 1 
- 1 


0.183 


+ 5 
- 5 




'^crit 


0.069 + 


1 
1 


0.051 


+ 1 
- 1 


0.179 


+ 5 
- 5 



Table 2: The decay constants (in lattice units) of the pseudoscalar and vector 
mesons. Also shown are the results for the combination fpy/Mp which in 
the heavy-quark limit is independent of the heavy-quark mass (up to mild 
logarithmic corrections) . 
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1 




Figure 1: The ratio of correlators defined in eq. (^) plotted versus t for 
Kh = 0.129, Ki = 0.14262. The curve represents the fit using timeslices 
15-22. 

We start the discussion of our results with the behaviour of the pseudoscalar decay constants 
as a function of the mass of the meson, with all dimensionful quantities given in lattice units. 
In the heavy-quark limit, the quantity fpy/Mp scales like 

Jp^/m^ = const. X [a,(Mp)]~^/'^° , Mp — > oo. (21) 

In order to detect possible deviations from this scaling law we plot in Fig. ^ the quantity^] 

l'(Mp) = («,(Mp)/«,(Mb))2/^°Z^Vp\/m^ (22) 

"'^The normalization factor as{MB)^^^^° is convenient when comparing these results with those obtained 
in the static theory. 
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as a function of 1/Mp. We approximate as{M) by 

27r 

where we take Aqcd = 200 MeV, and Pq = 11 — |nj, with n/ = in the quenched approx- 
imation. From the figure we see that <l>(Mp) increases as the mass of the heavy quark is 
increased (in agreement with the behaviour found using the Wilson action for the quarks 
B H) il)- order to quantify this behaviour, we fit <l>(Mp) to either a hnear or quadratic 
function of 1/Mp: 

or 

*(*W^c(l--^ + ^). (25) 



We have performed these fits twice; once using the values of fpy/Mp for all four values of k^, 
and once using those for only the smallest three k^s (i.e. for the heaviest three heavy-quark 
masses). The results of the fits are given in Table ^. We find that the non-scaling corrections 
are of 0(30%) for fp, and 0(10%) for fp, in agreement with previous results obtained using 
Wilson fermions 0,^3- From the quadratic fit to the data at all four heavy-meson masses 
we find, in physical units. 





= 0.45 


+ 2 

- 2 


+ 'I GeV^/^ 




= 0.84 


+ 11 

- 8 


+ ^3 GeV 




= 0.28 


+ 7 

- 9 


+ 'I GeV^. 



(26) 

The second error in eq. (p6D corresponds solely to the uncertainty in the scale. It should be 



mentioned that ignoring the residual logarithmic dependence of fpy/Mp on Mp makes the 
slope more pronounced. However it is clear from Fig. ^ and Table ^ that the logarithmic 
corrections to the scaling law can by no means account for the observed slope in fpy/Mp. 

We use the parameters of the fits in Table ^ to make our predictions for the values of the 
decay constants fn and fp- The results corresponding to the four fits are presented in 
Table |. From this table it is clear that there is a further systematic uncertainty in fp of 
about 11 MeV from extrapolating using either linear or quadratic fits. In contrast to this, 
since we interpolate to mp, the results for fp are very stable. It should be emphasised that 
choosing a different value for Aqcd (e.g. Aqcd = 250 MeV), or for the anomalous dimension 
(e.g., by taking nj = 4), changes the results by only about 1 MeV. 
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Linear Fit 


Quadratic Fit 




A 


B 


C 




E 


4 Kh's 
3 KhS 


0.089 + 2 
0.092 + 3 


0.199 + J 
0.216 + 


0.101 + 1 

~ 5 

0.102 + I 

— o 


0.31 t I 
0.33 + " 

— 5 


0.038 + J 
0.045 + 5^ 



Table 3: Values of the parameters of the linear and quadratic fits to the 
behaviour of the pseudoscalar decay constants with the mass of the mesons 
(as defined in the text). 





Linear Fit 


Quadratic Fit 




/d 


fB 






4 K/,'S 

3 K/i's 


185 + J ! ^5 

^gg + 4 + 41 


149 + ^ ! 52 
154 + 'I 

- 4 - 8 


185 !J 
185 +J 


160 ! g ! 
160 +! + 

- 7 - 7 



Table 4: Values of the decay constants Jb and in MeV, corresponding to 
the linear and quadratic fits. 



Taking the results from the quadratic fit using all four Kh values we find: 

fn = 185 + J + MeV (27) 

fB = 160 + ^ + II MeV (28) 

where we have included the uncertainty of 11 MeV from the extrapolations in the systematic 
error quoted for fs- We take the results presented in equations (|27| ) and ( p8D as our best 
estimates of the decay constants of the D and B mesons. 

In ref. [0] it was found useful to use the pion decay constant, /tt, to set the scale in the 
computations of the decay constants of heavy-light mesons. By calculating fo/ fw and fs/ fw 
it may be expected that some of the systematic errors cancel, since, in particular, the ratios 
are independent of Za- Our results for the decay constants obtained in this way are, as 
expected, close to the upper systematic error margins in Table ^. We find 

(^) X 132 MeV = 232 + MeV (29) 
(^) X 132 MeV = 201 + MeV. (30) 
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Figure 2: The data for $(Mp) plotted against the inverse meson mass. The 
open symbols denote points with ni < Kcrit, whereas full symbols denote 
those extrapolated to Kent- The solid line represents the linear fit to the 
chirally-extrapolated points using the three heaviest meson masses, whereas 
the dashed curve results from a quadratic fit to all four. 



Finally in this subsection, we present our results for and fs^- These are obtained by 
interpolating the measured values of the decay constants given in Table ^ to = k<j = 



0.1419 _ ^ |T^. The extrapolations in the heavy-quark masses are done as above. We find 



fn 

fB 



1.18 



1.22 



+ 2 

- 2 

+ 4 



(31) 
(32) 
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In physical units we obtain 



212 
194 



46 
7 

62 



MeV 
MeV. 



(33) 
(34) 



Recently the first measurement of f^^ has been made by the WA75 collaboration who 
found = (232 ± 45 ± 20 ± 48 )MeV. Our result is in good agreement with the measured 
value, and also with previous lattice calculations using Wilson fermions ^. 



2.2 Decay Constants of Vector Mesons 

In this subsection we present our results for the decay constants of heavy-light vector mesons. 
These are defined by 

m,\V)^e,-^ = Zymi:mV). (35) 
Jv 

where \ V) represents a state containing a vector meson V , with mass My, polarisation vector 
and decay constant fy- denotes the local lattice vector current, defined in eq. (||), 
with r = 7^, which has to be multiplied by the renormalisation constant Zy. The vector 
mass My is extracted from fits to the correlator 

c'y'y ^ EE(o|v}"(f,t)v;:^(o)|o) 

3=1 X 

exp(-MyLt/2)cosh(Mv(Lt/2-t)), (36) 



2My 

where is the jth spatial component of the smeared vector operator and Zys is defined 
through 

(0|l^^(0)|r) = e,Zys. (37) 

Fitting timeslices 14 < t < 23, symmetrized, we obtain the vector meson masses shown in 
Table |T]. In order to extract the matrix element of the local vector current we fit the ratio 

to a constant in the fitting interval 15 < t < 23. 

The results for 1/ fy for the twelve — i^i combinations are presented in Table |^, together 
with those obtained after extrapolation to the chiral limit. The chirally-extrapolated values 
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M 


Linear Fit 


Quadratic Fit 


oo 


1.02 + ] 

- 4 


1.09 ! \ 


(Mb + 3M^)/4 


0.93 + 3 


0.96 + ; 

— 5 


(Mb + 3M^)/4 


0.77 + I 


0.77 t 2 



Table 5: The quantity U{M) obtained from linear and quadratic fits. 



for fy Zy are now interpolated to the D* mass using a quadratic fit to the data at all four 
values of Kh, giving, 

/^3= 0.110 (39) 

This is slightly below, but still compatible with, earlier studies (e.g. 0) when the systematic 
error is taken into account. This result remains unaltered if a linear fit is used instead of a 
quadratic one. 



2.3 A Test of the Heavy Quark Symmetry 

In the heavy-quark limit, the decay constants of heavy-hght pseudoscalar and vector mesons 
are related by 



U(M)^fj^=(l + l^ + 0(l/M)y (40) 

where we take the heavy mass scale, M, to be the spin-averaged meson mass, M = {Mp + 

3My)/4. 

In order to test the predicted behaviour of U{M), we take the chirally-extrapolated values 
for both the pseudoscalar and vector decay constants, and fit 

t/(M).C/(M)/{l4^) (41) 

to either a linear or quadratic function of 1/M. The data together with the fits are shown 
in Fig. ^ and we display our results in Table ^ The perturbative values of Zj^ and Zy are 
used. 

The fact that [/(oo) is around one in Table | provides support for our parametrisations, in 
eqs. ( plf ) and (^Sf ), of the non-scaling behaviour of the decay constants for finite heavy-quark 
masses. 
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Figure 3: The quantity U (M) plotted against the inverse spin-averaged mass. 
Linear and quadratic fits are represented by the sohd and dashed curves, 
respectively. Also shown are the statistical errors of the extrapolation to the 
infinite mass limit. 



3 Decay Constants from the Simulation at = 6.0 

In this section we describe the results of a computation of the decay constants using the 
SW fermion action at (3 = 6.0 on a 16^ x 48 lattice. These results were obtained using 
36 configurations, with light-quark masses corresponding to k; = 0.1432,0.1440 and 0.1445. 
The corresponding light-light pseudoscalar and vector meson masses, and pseudoscalar decay 
constants, all in lattice units, are presented in Table ^. The values of the hopping parameter 
corresponding to the chiral limit and the strange quark mass are Kcrit = 0.14556 + I and 
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L 


TV 


p 


/ ^ A 

J TV / '-^ ri. 


0.1432 
0.1440 
0.1445 


0.386 + I 
0.311 + I 

— 5 

0.257 + I 


0.51 + \ 
0.47 + 2 
0.43 + 3 


0.088 + 3 
0.080 + I 
0.075 + ^ 

— 5 


/cerit = 0.14556 + I 




0.38 + J 


0.065 + I 

— D 



Table 6: Masses of light-light pseudoscalar and vector mesons, and the pseu- 
doscalar decay constants at /3 = 6.0. 

Kg = 0.1437 ~^ ^ respectively. Using the mass of the p meson to determine the value of 

the lattice spacing, we find = 2.0 + ^^eV, whilst using we find a"^ = 2.1 + JCeV. 
These two results are compatible, and below we will use the value 

a-^ = 2.0 + iGeV (42) 

to convert the results from lattice to physical units. 

We have computed the heavy-light correlation functions as series in (the hopping-parameter 
expansion [^), thus enabling us to obtain the decay constants at any value of the mass of 
the heavy quark, without explicitly computing the heavy-quark propagators. The decay con- 
stants are obtained by fitting to eq. ( pj]) and eq. (pOD, over the range 12 < t < 18 for both 
fits. We employ the Jacobi smearing algorithm with = 50, corresponding to a smearing 
radius of r = 4.2. 

In an attempt to improve our understanding of the discretisation errors, we have also com- 
puted the decay constants for the Wilson action at one value of the light-quark mass, using 
a subset of 16 of the 36 configurations. The comparison of the results for the two actions is 
presented in Subsection pl2| . 

3.1 Pseudoscalar Decay Constants 

In Fig. ^ we plot the chirally-extrapolated values of <l>(Mp) as a function of 1/Mp, for 11 
values of the heavy-quark mass. We fit the points corresponding to the five lightest meson 
masses (for which mq = 1/2(1/ Kh — l/z^crit) < 0.7, as was the case at /5 = 6.2) to eq. (pSf), 
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and this is shown as the sohd curve in the figure. For the coefficients of the fit we find: 



C = 0.18 + 3 ; D = 0.45 + ; E = 0.08 + g. (43) 
In physical units we obtain 

= 0.50 ! 9 ! GeV3/2 

Da~' = 0.91 ! 26 - '5 GeV 

Ea~^ = 0.32 + '° GeV' (44) 



in good agreement with the results at /5 = 6.2, quoted in eq. (|26|) . It should be noted that 
the inclusion of all 11 data points makes no significant difference to the fit. 

The values for the decay constants in physical units are: 

fD = 199 + ;j t Z MeV (45) 
fB = 176 MeV (46) 



Id 



24 - 15 

1.13 + I (47) 



^ = 1.17 + . (48) 

J B 

All these numbers are in good agreement with the corresponding results from the simulation 
at /5 = 6.2. Finally, for the ratios fo/ f-w and JbI J-k we obtain 

'^^^ X 132 MeV = 211 + MeV (49) 



^) X 132 MeV = 186 + MeV. (50) 



3.2 Comparison of Results Using Wilson and SW Actions 

For 16 of the configurations, we have computed the decay constants for the Wilson fermion 
action, again using the hopping-parameter expansion. We compute light-quark propagators 
at a single value of the hopping parameter, = 0.155, corresponding to a pseudoscalar- 
meson mass of 0.30 + J. This was chosen to match the SW pseudoscalar-meson mass of 

0.31 ~^ \ obtained at = 0.144 on the same set of configurations. 
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Figure 4: The chirally-extrapolated data for $(Mp) at /? = 6.0 plotted 
against the inverse meson mass. The sohd curve represents a quadratic fit 
to the points denoted by circles. Points represented by diamonds are not 
included in the fit. The dashed curves are the 68% confidence bounds on the 
fit. 



In Fig. I we plot fp^/lvl^ function of 1/Mp. Using the conventional normalisation of 
a/2k for the quark fields, we see a clear divergence between the results for the two actions for 
ntQ > 0.7; the Wilson results turn over and decrease. However, we note that uncorrelated 
fits of the Wilson points, at the lightest few meson masses, to eqs. (|^) and ( p5|) would 
yield coefficients of the 1/Mp term broadly consistent with previous Wilson analyses. 

The figure also shows the results obtained with the Wilson action, but using the normali- 
sation a/1 — Qk for the quark fields [^, |l^, where k = uqi^ and uq = l/(8Kcrit)- It has 
been suggested that this normalisation may absorb some of the discretisation errors |T^, 
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Figure 5: jp\/Mp for both the Wilson and SW actions. Diamonds denote 
points obtained with the Wilson action in the conventional normalisation, 
-\/2/t, whereas squares denote points normalised by Vl — 6k. Results using 
the SW action are represented by circles. The solid curves are quadratic fits 
in 1/Mp to fp\JMp for the Wilson action, with fields normalized by \/\ — Gk, 
and for the SW action. The dashed curve is to guide the eye. 



and indeed the corresponding results agree remarkably with those obtained using the SW 
action. This agreement provides considerable motivation for a theoretical study to investi- 
gate whether there is any formal connection between the ansatz above and the improvement 
programme initiated by Symanzik PT| . 
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4 /b in the Static Limit 



An alternative and complementary approach to heavy-quark physics using lattice QCD was 
proposed by Eichten [|1^]. This technique is based on an expansion of the heavy-quark 
propagator in inverse powers of the quark mass. In practice, one keeps just the leading term, 
given by (at zero velocity) 

SQ{x,t;6,0) = |^(t)e— «*i±^ + ^^(-t)e™«*^^| 6'^'Hx)V^{t,0), (51) 

where V^^t, 0) is the product of links from (0, t) to the origin, for example for t > 0, 

Pg(t, 0) = f/i(0, t - l)f/i(0, t - 2) ■ ■ ■ Ul{6, 0). (52) 
At sufficiently large times 

UMx,t)Aim ^ ^e-^'^\ (53) 



where is the improved axial current of eq. (10) with T = 7^75. Since the only dependence 
on rriQ in eq. (|53|) arises through the exponential factor in eq. (|5T|) , we deduce the scaling law 
that fp a/Mp is independent of the heavy-quark mass. Matching the result from the Heavy 
Quark Effective Theory with that in the full theory introduces the logarithmic corrections 
in eq. (0). The full scaling law is of the form 



/pv/m^ = const. [(a,(Mp))-2/A'(l + 0(a,)) + 0(1/Mp)] . (54) 

The objective of lattice computations is to determine the constant. We refer to the value of 
Jb obtained using eq. (0), but dropping the 0{1/Mb) corrections, as f^^^. 

We compute the two correlation functions, C^^ and C^^, defined byQ 

C'Ht) = Y.iO\AUx,t)Af (6,0)10) ^{Z'^fe-^''' (55) 

C'^s^t) = ^(0|A^(x,t)Af (0,0)|0) ^Z^ZV^^*, (56) 

X 

where AE is the (unphysical) difference between the mass of the meson and the bare mass 
of the heavy quark. The matrix element of the local operator is obtained from the two 
correlation functions C^'^ and C^^ as follows. By fitting C^^{t) to the functional form given 
in eq. (|55|) we obtain Zs (and AE). At sufficiently large times the ratio C^^{t)/C^^{t) —>■ 
Z^jZ^, so that Z^ can be determined. 



^In principle the behaviour of the correlation functions in eqs. (|5|) and (|5|) is given by a cosh (as in 
eq. (pj|)), however the contribution of the backward-propagating meson is negligible in the time intervals we 
will be considering. 
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Hi 


AE 






Z^ 




N = 140 


N = 110 


N = 140 


N = 110 


= 140 


N = 110 


N = 140 


0.14144 
0.14226 
0.14262 


0.59 + J 
0.57 + J 
0.56 + I 


141 + 

_ 

127 + ;^ 

119 


130 ! 
119 + 

112 + ' 


0.0125 + 3 
0.0121 + I 
0.0120 + 3 


0.0039 + J 
0.0038 + J 
0.0037 + J 


0.149 + I 
0.137 + I 
0.131 + ! 


0.142 + I 
0.130 + I 

— D 

0.125 + I 



Table 7: Values of AE, {Z^)^, Z^/Z^ and Z^ at the three value of ki. AE 
is obtained from the fit to C^^{t). 



4.1 Results at [3 = 6.2 

We now report on a computation of Z^**** at /3 = 6.2. The results presented here were obtained 
using a subset of 20 of the 60 configurations discussed in Section ^ at the three values of the 
light-quark mass. The values of Z^*'^* and were determined by extrapolating the results 
to Kcrit and Ks respectively. 

In view of the difficulty in isolating the ground state in correlation functions using the static 
effective theory, we have compared results obtained with different numbers of iterations of 
the Jacobi smearing algorithm |^ . For N less than about 80 the plateaus do not start until 



at least t = 7. In this paper we present our results obtained with = 110 and A^ = 140, 
corresponding to r = 5.9 and 6.4 respectively, where plateaus begin as early as t = 4 and 
hence statistical errors are smaller. 

In Fig. P(a) we show the effective masses obtained from C^^{t), and in Fig. ^(b) the ratio 
C^^{t)/C^^{t), both at Ki = 0.14226 and A^ = 140. Excellent plateaus are obtained, giving 
us confidence that the ground state has indeed been isolated. In Table |^ we present the 
results for AE, (Z^)"^, Z^/Z^ and Z^ at all three values of ki, from fits over the range 
5 < if: < 11, without symmetrization in Euclidean time. AE is obtained from the fit to 
C^'^{t) for A^ = 140; consistent values are obtained for A^ = 110. 

Extrapolating the results for Z^ in Table ^ to the chiral limit and to the mass of the strange 
quark we find: 

Z^ = 0.124 + ^ at Ki = Kcrit (57) 
Z^ = 0.140 at Ki = (58) 
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(a): ln(C^"(0/C^"(^+l)) 
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Figure 6: (a) The effective mass obtained from C^^{t), and (b) the ratio 
C^^{t)/C^^{t) at P = 6.2, Ki = 0.14226 and = 140. The sohd hnes 
represent fits from 5 < t < 11. 



when obtained using smeared interpolating operators with = 110 and 



0.117 



+ 7 



at Ki — Kcrit 



0.134 at Ki = Kg 



(59) 
(60) 



when using interpolating operators with = 140. 

When matching the static lattice theory to the full theory at a scale rrib, the factor required 
is [§§: 

as{a^^) r3 



7Stat 



Za\1 + 



37r 



- log a^ml - 2 



(61) 
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Za-, relating the axial current in the static lattice theory to the static continuum one, has 



been calculated in perturbation theory for the SW action [0, p!7[] : 

= 1 -0.127/ ~ 0.79. (62) 



The value of 0.79 on the right hand side of eq. ([6^ ) was estimated using the boosted coupling 
at /3 = 6.2. For the remaining factor in eq. (|6TD, we take = 5 GeV, given by eq. ( ]23|) 
with n/ = 0, and Aqcd = 200 MeV, yielding a number close to one (note that this is 
insensitive to small changes in nth). Thus Z^**** = 0.79 also. We find 

MeV (63) 

MeV, (64) 
when using the interpolating operators with = 110, and 

(7stat \ 
^) ^^^^ 

^ MeV (66) 
0.79 / 

when using those with = 140. The systematic errors quoted arise from the uncertainty in 




the scale. We take the results in equations ( p5| ) and (pB]) as our best values, and these give 
for the ratio: 

1.14 + 1 (67) 



ystat 

1b 



stat ■ — 3' 



4.2 Results at = 6.0 



We have performed a similar analysis on the 36 configurations at /3 = 6.0, using Jacobi 
smearing with = 50 and 150, corresponding to r = 4.2 and 6.2 respectively. The results 
obtained using the two smearing radii are consistent, and our best results are those at A^ = 50 
for which Zi^ = 0.211 ^ ^, yielding 



fT' = 286 : : ^ MeV (6^ 




/r = 323 : : ( ^ ) Mev. (69) 

xstat 

^ = 1-13 ! t- (70) 
Jb 
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(b): C^\t)/C^^{t) 








Figure 7: (a) The ratio C^^{t)/C^^ (t) and (b) C^^(t)/C^^(t) plotted 
against t for P = 6.0, = 0.1432, and = 50 iterations used in the 
Jacobi smearing algorithm. 



These results at j3 = 6.0 are consistent with those at /? = 6.2 presented in eqs. ([65|), ( ]66| ) 
and (|6l 



The results plotted in Fig. |(b) for the ratio C^^(t)/C^^(t) appear to have considerably 



smaller errors, and a clearer plateau, than in some recent studies |T0|, in spite of our limited 
statistics. We attribute this to the fact that we use the C^^ correlation function in which 
the smearing is performed at the source, rather than the C^^ correlation function in which 
the smearing is performed at the sinkQ. Of course, with sufficiently many configurations, 
the results are independent of this choice. However in the C^^ correlation function, the 



■^In both refs. g and it was in fact the C'^'^, and not the C^^ , correlation function which was 
computed 



25 



heavy-quark propagator is sampled at many spatial points, whereas in the C^^ correlation 
function only the heavy-quark propagator at = contributes. Thus it seems plausible that 
the statistical errors are considerably reduced using the C^^ correlation function. 

To check this hypothesis, we have computed the ratios C^^{t)/C^^(t) and C^^{t)/C^^(t) 
at /? = 6.0, with ni = 0.1432. The results are plotted in Fig. |^, and indeed confirm that 
there is an enormous improvement in precision when the correlation function C^^ is used. 
We believe that this, rather than the different method of smearing, is the reason for the 



relatively poor plateaus in ref. [10 



4.3 Discussion 

We begin with a comparison of the static and propagating results. Because we do not yet 
have static results for the full set of configurations at (3 = 6.2, we focus on a comparison at 
P = 6.0. In Fig. H we plot our results for the scaling quantity fpy/Mp{a{Mp)/a{MB))^^^^ 
from the simulation at (3 = 6.0 as a function of 1/Mp (in lattice units), together with 
our result for fp^^\/MB- The quadratic fit which we used to obtain our estimate for fp 
in Subsection |3.1| gives an intercept at 1/Mp = which is about 25% and two standard 
deviations below the static result; a similar discrepancy is observed at /3 = 6.2. There are 
a number of possible reasons for this, e.g. uncertainties in the renormalisation constants 
(which are different for the static and propagating quarks), residual discretisation errors in 
the simulation of the propagating quarks, and uncertainties in the various extrapolations. 

A better way of determining the consistency of the static and propagating results is to 
include the static result in the quadratic fit. Such a fit using the full correlation matrix at 
j3 = 6.0 gives a x^/dof of 1.5. This is still acceptable, and provides further evidence that 
using rotated operators with the SW action gives a sensible normalisation for propagating 
heavy-quark fields. From this fit we obtain fp = 220 ^ ^ 27 ^^V, which is 44 MeV higher 
than that obtained from the propagating points alone. 

In Table ^ we present the results for f]^^^ obtained by other groups, together with our values. 
Although at (3 = 6.0 the values of found by all groups and for both actions are in broad 
agreement, the different treatment of systematics leads to the spread of results in fp^^. It 



has been suggested that fp^^ decreases as a — > |jT2|. However, the agreement of the results 
obtained with the Wilson and SW actions at (3 = 6.0, together with consistency between 
our results at /? = 6.0 and P = 6.2, suggests that the discretisation errors are small. 
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Ref. 


Action 




a^i [GeV] 


vstat 


[MeV] 




11 




Wilson 


5.9 


1.75 


0.79 


319 ± 11 




i 




Wilson 


6.0 


2.0 ±0.2 


0.8 


310 ±25 ±50 




i 




Wilson 


6.0 


2.2 ±0.2 


0.8 


366 ± 22 ± 55 




10 


1 


Wilson 


6.0 


0"ii 1 r\ r\r' i r\i/~i 

2.11 ± 0.05 ± 0.10 


0.8 


350 ± 40 ± 30 




10 


1 


SW 


6.0 


2.05 ± 0.06 


0.89 


370 ± 40 


This Work 


SW 


6.0 


2.0 


0.78 


„„„ +8+67 

286 ! + 2 


This Work 


SW 


6.2 


2.7 


0.79 


rjp-O + 16 + 105 
^^'^ - 15 - 14 




i 




Wilson 


6.3 


3.21 ±.09 ±.17 


0.69 


235 ±20 ±21 




12 


1 


Wilson 


5.74, 6.0, 6.26 


1.12, 1.88, 2.78 


0.71(8) 


230 ± 22 ± 26 



Table 8: Compilation of lattice results for f\ 



5 Conclusions 

In this paper we have carried out an extensive study of the decay constants of heavy-light 
mesons using the SW action for the quarks. The use of the SW action confirms the large, 
negative OiXjMp) corrections to the scaling law fpy/Mp ~ constant at the mass of the 
charm quark and the significant corrections at the mass of the b quark, previously observed 
with the Wilson action. However, from our comparison of results for the Wilson and SW 
actions at /5 = 6.0, we observe clear evidence that the ^/2k normalisation of the Wilson 
quark fields fails for large quark mass. This failure is presumably due to large 0{mQa) 
effects. It has been suggested that such effects may largely be absorbed by the use of a 
different normalization We find that such a normalization yields results in agreement 
with those obtained using the SW action with rotated operators and the \/2k normalization 
for the quark fields. 

Our best estimates of fo and fs are 

= 185 + J (stat) + (syst) MeV (71) 

fs = 160 + J + II MeV, (72) 

obtained using propagating quarks at /5 = 6.2. Our analysis at /? = 6.0 yields entirely 
consistent results. The latter analysis also suggests that including the static result in the 
fits is likely to increase the value of Jb by around 40 MeV. 



27 




1 2 

1/Mp 



Figure 8: Za^{Mp) at /? = 6.0 from the simulation with propagating quarks 
(open symbols) and the static theory (cross). The dashed curve is the fit to 
the open circles, with the parameters of eq. (|43|) ; the square is the intercept 
at 1/Mp = 0. The solid curve is the fit with the static point included. 

The most urgent extension of this work is to determine the B parameter for B^-B^ mixing, 
since it is the combination fBy/Bp which is directly relevant for phenomenological studies 
of the mixing and of CP- violation. A recent simulation with Wilson fermions found Bp = 
1.16 ± 0.07 0, and it is important to confirm this result with the improved action. 
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